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QCD at high density/temperature 
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I summarize recent work on QCD at high temperature and density. 



1. INTRODUCTION 

QCD now occupies an established position 
as the well-attested theory of strong-interaction 
physics. Other contributions to these proceedings 
(those of S. Frixione, L. Lellouch, and K. Long) 
survey the agreement between evidence from par- 
ticle collider experiments and perturbative and 
lattice QCD calculations. 

However, there is more to QCD than collisions 
between pairs of particles. QCD should also be 
able to tell us about the thermodynamics of mat- 
ter in the realm of extraordinarily high temper- 
atures (> 100 MeV) and densities at which it 
comes to dominate the behavior. This regime 
has real physical importance: the whole universe 
was hotter than 100 MeV for the first crucial 
microseconds of its history, and in the current 
epoch our galaxy has a large population of neu- 
tron stars consisting of matter squeezed beyond 
nuclear density by gravitational forces. However, 
only in the last few years have these regimes be- 
gun to be probed experimentally in heavy ion col- 
lisions and astrophysical observations of neutron 
stars, and we are still working towards achieving 
a proper theoretical understanding of them. 

In this paper I will survey the progress that 
has been made in the last few years in qualita- 
tive predictions and quantitative calculations of 
the properties of QCD at very high temperature 
and net baryon density. I will not discuss theory 
that is specifically tailored to heavy ion collisions, 
but refer the reader to existing review articles . 
Other more tightly focused reviews of QCD at 
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Figure 1. A conjectured phase diagram for QCD 
in the real world. The CFL region may alterna- 
tively extend down to the nuclear phase, with no 
intermediate phase. Depending on the strength 
of instanton interactions, the CFL phase may in- 
clude K'^ condensation. 



high density [|| or high temperature |^ may also 
prove useful for details that go beyond the scope 
of this summary. 

2. Phases of QCD 

Fig. shows a conjectured phase diagram for 
QCD. Along the horizontal axis the temperature 
is zero, and the density rises from the onset of 
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nuclear matter through the transition to quark 
matter. Compact stars are in this region of the 
phase diagram, although it is not known whether 
their cores are dense enough to reach the quark 
matter phase. Along the vertical axis the temper- 
ature rises, taking us through the crossover from 
a hadronic gas to the quark gluon plasma. This 
is the regime explored by high-energy heavy-ion 
colliders such as RHIC (see the contribution by 
J. Harris to these proceedings). At high tempera- 
ture and moderate density we find the first-order 
chiral phase transition ending at the chiral critical 
point which could perhaps be detected [5[| in 
medium-energy heavy ion colliders (see SectT^. 

For some parts of the phase diagram reliable 
calculations can be performed. Standard pertur- 
bative methods can be used in the outer regions of 
high density (/i > 10® MeV |^) and temperature 
(T > 10^®'^ GeV 0), and resummation methods 



can extend this down to T > 3Tc (Section 4.2). 
At zero density and T < 4Tc, lattice methods give 
good results, and ways have recently been found 
to extend lattice methods to fi < T. However, 
this leaves much of the most interesting phase 
structure in regions where controlled calculations 
are not yet possible. In the following sections, I 
will describe the progress that has been made in 
understanding QCD at high density, in achieving 
first-principles calculations at high temperature, 
in lattice calculations at high temperature and 
low density, and in gaining control of field theo- 
ries out of thermal equilibrium. 



all, QCD binds quarks together to form baryons) 
and it was shown by Bardeen, Cooper, and Schri- 
effer (BCS) @ that if there is any channel in 
which the interaction is attractive, then there is 
a state of lower free energy than a simple Fermi 
surface. That state arises from a complicated 
coherent superposition of pairs of particles (and 
holes) — "Cooper pairs". 

Heuristically, the BCS argument is this. The 
Helmholtz free energy is F = E — ^N, where E 
is the total energy of the system, /i is the chem- 
ical potential, and N is the number of fermions. 
The Fermi surface is defined by a Fermi energy 
Ep — fi, at which the free energy is minimized, 
so adding or subtracting a single particle costs 
zero free energy. Now switch on a weak attrac- 
tive interaction. It costs no free energy to add 
a pair of particles (or holes), and the attractive 
interaction between them then lowers the free en- 
ergy of the system. Many such pairs will therefore 
be created in the modes near the Fermi surface, 
and these pairs, being bosonic, will form a con- 
densate. The ground state will be a superposition 
of states with all numbers of pairs, breaking the 
fermion number symmetry. 

Since pairs of quarks cannot be color singlets, 
the resulting condensate will break the local color 
symmetry SU{3)co\ot- We call this "color super- 
conductivity" . Note that the quark pairs play the 
same role here as the Higgs particle does in the 
standard model: the color-superconducting phase 
can be thought of as the Higgs phase of QCD. 



3. High density QCD 

It is difficult to make accurate predictions for 
high density QCD because lattice gauge theory is 
stymied by the complexity of the fermion deter- 
minant (see Sect. H). However, we can borrow 
some intuition from condensed matter physics, 
which tells us to expect the very interesting phe- 
nomenon of color superconductivity [^,^ QCD 
is asymptotically free — the interaction becomes 
weaker as the momentum grows — so at suffi- 
ciently high density and low temperature, there 
is a Fermi surface of almost free quarks. The in- 
teractions between quarks near the Fermi surface 
are certainly attractive in some channels (after 



3.1. Calculation methods 

To decide whether or not fermions condense in 
the ground state, one can explicitly construct a 
wavefunctional with the appropriate kind of pair- 
ing, and use a many-body variational approach 
to find the pairing strength by minimizing the 
free energy. Alternatively but equivalently, one 
can formulate the problem field-theoretically, in 
terms of mean-field Schwinger-Dyson equations 
for a quark self energy that includes "anomalous" 
(qq) terms as well as the usual {qq) terms. In ei- 
ther case the problem reduces to solving a self- 
consistency equation, the gap equation, to find 
the self energy. If the only solution is zero, there 
is no condensation in that channel. If not, there 
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can be condensation, but it may just be a local 
minimum of the free energy. There may be other 
solutions to the gap equation, and the one with 
the lowest free energy is the true ground state. 

There are several possible choices for the in- 
teraction to be used in the Lagrangian from 
which we derive the gap equation. At asymp- 
totically high densities {fi > 10® MeV ||) QCD 
is weakly coupled, so one gluon exchange is ap- 
propriate. Such calculations |11,1^ demonstrate 
from first principles that color superconductivity 
occurs in QCD. However, the density regime of 
physical interest for neutron stars or heavy ion 
collisions is up to a few times nuclear density 
(/X < 500 MeV) , so one must use a phenomenolog- 
ical interaction that can be argued to capture the 
essential physics of QCD. The usual choice is an 
NJL model, with a four-fermion interaction based 
on known QCD interactions (gluon exchange, or 
instantons) whose strength can be normalized 
to reproduce known low-density physics such as 
the chiral condensate, and then extrapolated to 
the desired chemical potential. Such calculations 
|8||^, [l3| , |l4| find that the color superconducting 
gap is insensitive to the form of the interaction, 
and rises to values around 100 MeV. Various 
other calculations that include competition be- 
tween chiral symmetry breaking and quark pair- 
ing confirm this . Very recently, the first lat- 
tice calculation of a 3-1-1 dimensional NJL model 
with chemical potential has been performed, and 
there are indications of diquark condensation [|l6j . 

The favored pairing pattern at high densities, 
where many strange quarks are present, is "color- 
flavor locking" (CFL), which will be described be- 
low. At lower densities there may be a direct tran- 
sition from nuclear matter to CFL quark matter 
(see Section \i.6^ . Alternatively, at intermediate 
densities there may be another phase. One pos- 
sibility is the "2SC" phase which has u-d pair- 
ing only, but this now seems unlikely to occur in 
equilibrated quark matter (Section |3.3| ). Other 
poss ibil ities include crystalline pairing (see Sec- 
tion 



3.7), or some sort of single-flavor pairing |lj 



3.2. Three flavors: Color- flavor locking 
(CFL) 

In QCD with three flavors of massless quarks 
the Cooper pairs cannot be flavor singlets, and 
both color and flavor symmetries are necessar- 
ily broken |l|]. Both NJL H,|l8|jl9| and gluon- 
mediated interaction calculations |12| agree that 
the attractive channel exhibits a pattern called 
color-flavor locking (CFL) |^, 



(1) 



[5[/(3)coior] X SU{3)l X SU{3)rxU{1)i 

3 [C/(1)q] 
^ SU{3)c+l+r xZ2 

' 3 [C/(1)q] ' 

(color indices a, (3 and flavor indices i, j run from 
1 to 3). The Kronecker deltas connect color in- 
dices with flavor indices, so that the condensate is 
not invariant under color rotations, nor under fla- 
vor rotations, but only under simultaneous, equal 
and opposite, color and flavor rotations. Since 
color is only a vector symmetry, this condensate 
is only invariant under vector flavor -f color rota- 
tions, and breaks chiral symmetry. The effect of 
including a sufficiently large strange quark mass 
is to unpair the strange quark pl| , p^ , p2[ , and 
at lower values it may induce a flavor rotation 
of the condensate kno wn a s "kaon condensation" 
(Ref. H and Section |3Tll) . 



The features of this pattern of condensation are 

• The color gauge group is completely broken. 
All eight gluons become massive. This en- 
sures that there arc no infrared divergences 
associated with gluon propagators. 

• All the quark modes are gapped. The nine 
quasiquarks (three colors times three fla- 
vors) fall into an 8 © 1 of the unbroken 
global SU{3), so there are two gap parame- 
ters. The singlet has a larger gap than the 
octet. 

• A rotated electromagnetism ( "Q" ) survives 
unbroken. It is a combination of the original 
photon and one of the gluons. 

^Interestingly, such a pattern was at one time considered 
as a possibility for symmetry breaking at zero densitybu] . 
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• Two global symmetries are broken, the chi- 
ral symmetry and baryon number, so there 
are two gauge-invariant order parameters 
that distinguish the CFL phase from the 
QGP, and corresponding Goldstone bosons 
which are long-wavelength disturbances of 
the order parameter. If a quark mass were 
introduced then it would explicitly break 
the chiral symmetry and give a mass to the 
chiral Goldstone octet, but the CFL phase 
would still be a superfluid, distinguished by 
its baryon number breaking. 

• The symmetries of the 3-flavor CFL phase 
are the same as those one might expect for 
3-flavor hypernuclear matter so it is 
possible that there is no phase transition 
between them. 

3.3. Absence of 2SC 

Until recently it was thought that if there were 
not enough strange quarks present in quark mat- 
ter, as might happen at intermediate density, the 
favored pairing pattern would be 2SC, in which 
the u and d quarks of two colors pair, leaving 
the s quarks and third color quarks unpaired. 
But in a real neutron star we must require elec- 
tromagnetic and color neutrality [ p4|j2^ (ignor- 
ing charge-separated phases, see Section 3.6). It 
turns out that this penalizes the 2SC phase rel- 
ative to the CFL phase ||2^. The reason is that 
the CFL phase has already paid most of the cost 
of neutrality, since it brings the u, d, and s Fermi 
surfaces close together p2| , p5|p^ , and the 2SC 
phase only pairs 4 of the 9 quark colors and fla- 
vors, so it has much less pair binding energy than 
the CFL phase. 

The arguments made in Ref. are model- 
independent, based on simplified assumptions 
about the dependence of the constituent strange 
quark mass Ms on /i and expanding the free en- 
ergy in powers of Mg/fi. The NJL calculation of 
Ref. handles Ms ^ /i and includes the cou- 
pling between the chiral condensate and quark 
condensate gap equations. The net result is that 
once neutrality is imposed, there is no (or very lit- 
tle [|^) density range in which 2SC is the phase 
with the lowest free energy. 



3.4. Signatures in compact stars 

The only place in the universe where we expect 
sufficiently high densities and low temperatures 
is compact stars, also known as "neutron stars" , 
since it is often assumed that they are made pri- 
marily of neutrons. A compact star is produced in 
a supernova. Typical compact stars have masses 
close to 1.4M0, and are believed to have radii 
of order 10 km. although uncertainty about the 
equation of state leaves us unsure of the radius 
and the maximum density attained in the core. 
During the supernova, the core collapses, and its 
gravitational energy heats it to temperatures of 
order 10^^ K (tens of MeV), but it cools rapidly 
by neutrino emission. Within a few minutes its 
internal temperature T drops to 10^ K (100 keV). 
Neutrino cooling continues to dominate for the 
first million years of the life of the star. (For de- 
tailed treatments and and reviews see Ref. pTf). 

Color superconductivity affects the equation of 
state at order (A/^)^. It also gives mass to ex- 
citations around the ground state: it opens up a 
gap at the quark Fermi surface, and makes the 
gluons massive. One would therefore expect its 
main physical manifestations to relate to trans- 
port properties, such as mean free paths, conduc- 
tivities and viscosities. In the following sections 
we will survey some of these manifestations. 

3.5. Mass-radius relationship 

Although the effects of color superconductivity 
on the quark matter equation of state are sub- 
dominant, they may have a large effect on the 
mass-radius relationship. The reason for this is 
that the pressure of quark matter relative to the 
hadronic vacuum contains a constant (the "bag 
constant" B) that represents the cost of disman- 
tling the chirally broken and confining hadronic 
vacuum. 



(2) 



If the bag constant is large enough so that nuclear 
matter is favored (or almost favored) over quark 
matter at /i ^ 320 MeV, then the bag constant 
and /i^ terms almost cancel, and the supercon- 
ducting gap A may have a large effect on the 
equation of state and hence on the mass-radius 
relationship of a compact star p8[. 
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3.6. Interfaces and mixed phases 

If color superconducting quark matter exists in 
compact stars, we expect it to occupy a core re- 
gion, surrounded by a mantle of nuclear matter. 
There may be a sharp interface between the two, 
but it is well known that a mixed phase region 
is also possible , given that there are multiple 
chemical potentials in nuclear matter: for baryon 
number, electric charge, and also color. A mixed 
phase region would have distinctive characteris- 
tics. Its transport properties are very different 
from those of the uniform CFL state: neutrino 
mean free paths, which are hundreds of meters in 
the CFL phase at temperatures of a few MeVpO[, 
are short in the mixed phase due to coherent scat- 
tering 

The nature of the boundary between nuclear 
and CFL quark matter has been studied p2[ , 
and it was found that a mixed phase only oc- 
curs if the surface tension of the interface is less 
than about 40 MeV/fm^ = 0.2 x (200 MeV)^, 
a fairly small value compared to the relevant 
scales Aqcd ~ 200 MeV, /i - 400 MeV. This 
means that there may well be a sharp interface 
rather than a mixed phase, with a baryon and 
energy density discontinuity of about a factor of 
two across it. Such an interface may modify the 
mass vs. radius relationship for compact stars 
with quark matter cores. It may also have qual- 
itative effects on the gravitational wave profile 
emitted during the in-spiral and merger of two 
compact stars of this type. Finally, it will affect 



the r-mode spectrum (see Section 3.10), and the 
damping forces to which r- modes are subject. 

3.7. Glitches and the crystalline color su- 
perconductor 

There is a crystalline form of quark pairing (the 
"LOFF" phase) that occurs when two differ- 
ent types of quark have sufficiently different Fermi 
momenta (because their masses or chemical po- 
tentials are different) that BCS pairing cannot 
occur 110 

Such situations are likely to be generic in na- 
ture, where, because of the strange quark mass. 



combined with requirements of weak equilibrium 
and charge neutrality, all three flavors of quark in 
general have different chemical potentials. 

The phenomenology of the crystalline phase 
has not yet been worked out, but recent calcu- 
lations using Landau-Ginzburg effective theory 
indicate that the favored phase may be a face- 
centered cubic crystal with a reasonably 
large binding energy. This raises the interesting 
possibility of glitches in quark matter stars. 

Glitches are sudden jumps in rotation fre- 
quency O of a pulsar. They are thought to arise 
from pinning of superfluid vortices in the nuclear 
matter core to the lattice of nuclei that forms the 
crust As the pulsar spins down, the vor- 

tices move outwards, and eventually overcome the 
pinning forces, and rearrange themselves. Un- 
til recently it was thought that quark matter 
stars would not glitch, because quark matter was 
thought to have no crust or other rigid structure 
that could pin the vortices. Now, however, we 
can envisage an alternative scenario for glitches 
(supported only be dimensional estimates at this 
point) where a superfluid CFL core is surrounded 
by a crystalline layer of quark matter in a LOFF 
phase. 

3.8. Cooling by neutrino emission 

For its first million or so years, a compact star 
cools by neutrino emission. There are many un- 
certainties in determining the temperature of a 
compact star (from its X-ray spectrum) and its 
age (inferred from the spindown rate) but ob- 
servers are now able to map out cooling curves 
||3|,^,Q. The cooling rate is determined by 
the heat capacity and emissivity, both of which 
are sensitive to the spectrum of low-energy ex- 
citations, and hence to color superconductivity 

il,|i|,il,io|. 

In the CFL phase, all quarks and gluons have 



^ At large number of colors there is the possibility of a sep- 
arate "chiral crystal" phase [ p^ , which is not a superfluid, 
so it will not lead to glitches. 



gaps A ^ T, electrons are absent [25 , and the 
pseudo-Goldstone bosons have masses of order 
tens of MeV ji^ , so the transport properties are 
dominated by the only true Goldstone excitation, 
the superfluid mode arising from the breaking of 
the exact baryon number symmetry. This means 
that CFL quark matter has a much smaller neu- 
trino emissivity and heat capacity than nuclear 
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matter, and hence the coohng of a compact star 
is hkely to be dominated by the nuclear mantle 
rather than the CFL core (^[||,|||. A CFL core 
is therefore not expected to be detectable by cool- 
ing measurements, although it has been suggested 
that it might slow the first « 100 years of cooling 

il- 

Other phases such as 2SC or LOFF give large 
gaps to only some of the quarks, leaving the oth- 
ers with no gaps (or very small gaps in weakly 
attractive channels) so they may allow fast cool- 
ing, which is severely constrained by the latest 
space-based X-ray observations . The cooling 
would proceed quickly, then slow down suddenly 
when the temperature fell below the smallest of 
the small weak-channel gaps. If no sign of such 
behavior is seen as our observations of compact 
star temperatures improve then we will have to 
conclude that either these phases do not occur, 
or their gaps are much larger than we expected. 

3.9. The neutrino pulse at birth 

In the first minute of a supernova, the inner 
regions cool rapidly, radiating off much of their 
energy as neutrinos, whose detailed spectrum as 
a function of time is determined by the neutrino 
diffusion properties of the protoneutron star, and 
can be detected in terrestrial experiments such as 
SuperKamiokande . 

Initial calculations by Carter and Reddy |Q 
found a striking signature: the core released few 
neutrinos when it was close to the critical temper- 
ature for color superconductivity T w T^, then a 
late blast of neutrinos as their mean free path 
grew larger than the size of the star at T ^ Tc. 
This picture is complicated by scattering of neu- 
trinos by Goldstone bosons Also, the tem- 
poral variation of neutrino emission from the core 
over the first few seconds is liable to be blurred 
by interactions with the surrounding material. 
Detailed numerical simulations of core collapse, 
including the effects of quark pairing, are now 
needed. 

3.10. r-mode instabiUty 

An "r-mode" (short for "rotational mode") 
is a bulk flow in a rotating star that ra- 
diates away energy and angular momentum in 



the form of gravitational waves (for a review, 
see Ref. @). If the rotation frequency of the 
star is above a critical value, the system becomes 
unstable to r-modes and quickly spins down to 
the critical value, at which point the r-modes 
are damped out. The critical frequency depends 
on the sources of damping, such as shear and 
bulk viscosities, and also friction at the inter- 
face between the r-mode region and adjacent re- 
gion. One can therefore rule out certain models 
for compact stars on the grounds that they have 
such low damping that they could not support the 
high rotation rates observed in pulsars: r-mode 
spindown would have slowed them down. 

It is generally believed that ordinary nuclear 
matter is not ruled out in this way [Q , but Mad- 
sen has shown that color superconductivity, 
by creating gaps in the quark excitation spec- 
trum, suppresses the viscosities by factors of order 
exp(— A/r), encouraging r-mode spindown. Fur- 
thermore, since CFL matter contains no electrons 
p5|] , there is no electrostatic cushion to support 
a crust of normal matter, which would otherwise 
have been an additional source of friction. Mad- 
sen found that for a compact star made entirely 
of quark matter in the CFL phase, even a gap as 
small as A = 1 MeV is ruled out by observations 
of millisecond pulsars. 

It remains to extend this calculation to the 
more generic picture of a quark matter core sur- 
rounded by a nuclear mantle. The friction at the 
core-mantle interface may be enough |^,^ to 
quash the r-modes. Finally, quark matter may 
contain a region of LOFF crystal, and the r- 
modes could be damped at the edges of that re- 
gion rather than at the crust. 

3.11. Other phenomenology 

Many other physical phenomena associated 
with color superconductivity in quark matter 
have been investigated. Small chunks of quark 
matter ( "strangelets" ) may exist, and standard 
calculations must be modified in that context to 
account properly for finite volume effects [S^ . 
Strangelets in the CFL phase will have a differ- 
ent mass-charge relationship from those made of 
unpaired quark matter . At high enough den- 
sities U{1)a breaking is suppressed, leading to do- 
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main walls bounded by strings |p2| . 

There are many interesting questions about the 
influence of the strange quark mass on color su- 
perconductivity. When rUs = the favored phase 
is CFL, as discussed in Section 3^ . As the strange 
quark mass is increased from zero the CFL pair- 
ing becomes harder and harder to maintain, and 
at some point it ceases to be the ground state 
21 1 . But before that the system may respond to 
the finite by rotating the quark-quark con- 
densate in flavor space, in a direction that re- 
duces the number of strange quarks. This can be 
understood as condensation of kaons 1 23 , where 
by "kaon" we mean the pseudo Goldstone bosons 
that are associated with the breaking of the chi- 
ral flavor symmetry by the CFL condensate. The 
resultant CFL-K° phase would break isospin and 
may also include charged kaons in the first sec- 
onds after the supernova, before neutrinos escape 
the star [|3[|5|. 

It is unclear whether kaon condensation will oc- 
cur in dense quark matter cores of compact stars. 
Instanton and electromagnetic interactions tend 
to disfavor CFL-K" over ordinary CFL, and the 
magnitude of instanton effects at the densities of 
interest is very uncertain But CFL-K° has 
many interesting features. It supports topological 
or near-topological solitons: superconducting K- 
strings and domain walls |^ whose possible 
observable manifestations are yet to be investi- 
gated. 



4. High temperature QCD 

We now turn from the horizontal axis of the 
QCD phase diagram (Fig. |l|) to the vertical axis, 
and study QCD at high temperature but low 
(net) baryon number density. The standard tech- 
niques are lattice gauge theory, perturbation the- 
ory, and a hybrid of the two that is generally 
known as "dimensional reduction" . 

4.1. Lattice gauge theory 

Lattice gauge theory is a brute-force evalua- 
tion of the functional integral over gauge and 
quark field configurations that is the theoretical 
definition of QCD. It is therefore a fully non- 
perturbative first-principles calculation of the 



predictions of QCD. To be numerically tractable, 
the integral must be formulated in Euclidean 
time, 

Z = J VAV^ViP exp{-^M[A]iP) 

exp(-S'gauge[^]) (3) 

= /23Adet(M[A])exp(-S'gaugc[A]) 

where the fermionic matrix AI is the usual 
fermion action 



M[A] 



m 



(4) 



In general, lattice gauge calculations require large 
computational resources. A typical calculation is 
formulated in a box of size F « (3 fm)'^, with 
lattice spacing a ~ 0.1 fm, and with up and down 
quark masses rriu^d ^ To obtain physi- 

cally relevant results it is necessary to extrapo- 
late a ^ 0, V ^ oo, rriu^d Wg. In addition, 
not all quantities are accessible to lattice calcu- 
lations: spectral functions and transport coeffi- 
cients are defined in Minkowski space, and can 
only be obtained by some sort of analytic contin- 
uation from the lattice results |5^. In spite of 
these challenges, lattice gauge calculations have 
been very successful in obtaining predictions for 
the equation of state of QCD at high tempera- 
ture. 

In a lattice calculation, thermal averages are 
obtained by setting the length of the Euclidean 
time direction to the inverse temperature. To 
probe physics near the deconfinement tempera- 
ture Tc ~ 200 MeV ~ (1 fm)"-"^ therefore requires 
a lattice spacing a < 0.2 fm. To probe much 
higher temperatures would require a correspond- 
ingly smaller lattice spacing, increasing the nu- 
merical effort. Thus calculations are so far lim- 
ited to T < 4Tc. (This restriction may be eased in 
the future by the introduction of anisotropic lat- 
tices, which are finer in the ^-direction.) State-of- 
the-art lattice calculations use improved actions, 
which include extra terms to compensate for the 
0(a^) discretization errors, and yield continuum- 
like results from relatively coarse (a < 0.2 fm) 
lattices. This is crucial because the cost of a lat- 
tice calculation rises as a very high power of l /g . 
Using improved staggered |^ and Wilson |60| 
actions the equation of state of QCD has been 
calculated and is shown in Fig. 0. 
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Figure 2. The pressure in QCD with various num- 
bers of dynamical quarks, as a function of tem- 
perature. Arrows show the corresponding nonin- 
teracting quark-t-gluon gas pressure psB- For the 
(2-|-l)-flavor calculation the light quark mass was 



The fact that at T ~ AT^ the pressure is about 
15% less than the pressure psB of free quarks and 
gluons shows that there are still important inter- 
actions at that temperature. The major challenge 
for other approaches is to confirm and explain 
this. 

4.2. Perturbation theory 

Perturbation theory is generally applicable 
when the coupling is small. One might expect 
that quarks and gluons at temperature T would 
be described by QCD with running coupling g{T), 
which becomes small at high temperature, caus- 
ing perturbation theory to converge quickly. In 
fact, this does not in general happen. For exam- 
ple, attempts to calculate the free energy to order 
0{g^T'^) 1^ 61,6|] find that there are large cor- 
rections with alternating signs, making the naive 
perturbative series useless for physically interest- 



ing temperatures (for illustration, the g"^ predic- 
tion is shown in Fig. ^). As we will see below, 
there are resummations of the perturbative series 
that greatly improve its convergence. There are 
also now results to leading order in g for transport 



coefficients, whose expansion contains powers of 
log(.g) and therefore converges even more slowly. 

4.2.1. Equation of state 

There have been various proposals for reorga- 
nizing the high-temperature perturbative series 
for the pressure, in order to improve its conver- 
gence. They involve resumming the hard ther- 
mal loops (HTL), typically via the Cornwall- 
Touboulis-Jackiw two-particle irreducible (2PI) 
effective action Q , in which the complete prop- 
agator is used as an infinite set of variational pa- 
rameters |6j,^. This has been used in approxi- 
mately self-consistent HTL resummation calcula- 
tions of the entropy by Blaizot, lancu and Rebhan 
(BIR) and of the pressure by Peshier |^ . In 
a different approach, Andersen, Braaten, Petitgi- 
rard, Strickland (ABPS) used the HTL effective 
action with the Debye screening mass as a vari- 
ational parameter The results of such cal- 
culations are shown in Fig. |^. The dashed lines 
show 0(5^) naive perturbation theory, to illus- 
trate the poor convergence of the perturbative se- 
ries. (The 0{g^) and @{g'^) predictions lie below 
the free gas result.) The theoretical uncertainty 
comes from varying the renormalization scale by 
a factor of 2 around the first Matsubara frequency 
2ttT. The grey bands show the results of lattice 
calculations . Finally, the range of predictions 
of the BIR and ABPS hard thermal loop methods 
are shown. The BIR prediction has a bigger theo- 
retical uncertainty because they choose a scale to 
separate hard and soft contributions, and vary- 
ing that scale contributes to the error. BIR are 
currently working to reduce this uncertainty. 

We conclude that there are (perhaps some- 
what ad hoc) resummations of perturbation the- 
ory that indicate that the pressure should slowly 
approach the free QGP value Psb from below. 
This is in agreement with, but not yet as accu- 
rate as, the lattice calculations. At this point 
it appears difficult to push these calculations to 
higher order and thereby judge their validity. 

4.2.2. Transport properties 

The photon emission rate from a quark-gluon 
plasma has been perturbatively estimated by var- 
ious groups ||7l|,[72|j7^ to leading order in the 
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Figure 3. Results [Ty] of calculations of 

the pressure of hot gluons, from lattice calcu- 
lations (Grey bands) and &{g^) naive perturba- 
tion theory (dashed lines). The dot-dashed lines 
show the range of predictions from the HTL re- 
orderings of the perturbative series formulated by 
Blaizot, lancu, and Rebhan (BIR), and Andersen, 
Braaten, Petitgirard, and Strickland (ABPS). 



electromagnetic and strong coupling constants. 
However, there are complications due to lead- 
ing logs that arise from multiple scattering of 
co-linear bremsstrahlung (Landau-Pomeranchuk- 
Migdal (LPM) effect). Only recently have these 
been fully taken into account by Arnold, Moore, 
and Yaffe @. They found that the soft (fc < 2T) 
photon emission rate behaves as dT/dk oc fc^^/^, 
which is less infrared singular than the result 
dT/dk oc k~^ obtained in earlier calculations that 
neglected the LPM effect [73,73]. Similar meth- 



ods are applicable to the dilepton rate, which is 
an important observable in heavy-ion collisions, 
and also to transport properties such as the shear 
viscosity, electrical conductivity, and flavor diffu- 
sion constant, which are in principle important 
input parameters to phenomenological models of 
the early stages of a heavy-ion collision. We can 
therefore look forward perturbative calculations 
of these quantities to leading order in g in the 
near future. 



4.3. Dimensional reduction 

The intuition behind dimensional reduction is 
related to the poor behavior of the perturbative 
expansion in QCD at high temperature, which 
can be understood by considering the T —> oo 
limit of the pure gauge theory. At non-zero T 
the Euclidean time size is finite, so the glue fields 
can be Fourier-analyzed into Matsubara modes, 
with frequencies w = 0, 2ttT, AttT, etc. As T ^ 
oo, the Euclidean time size goes to zero, and the 
uj = Q Matsubara "static mode" dominates. The 
theory reduces to a three-dimensional theory of 
gluons (the spatial components Ai) and a scalar 
that transforms as a color octet {Aq) p^jpG 77 



C = ^tiF^.+tiiD„Aof+mlti-Al + XA{trAlf{5) 



This effective theory is confining, which indi- 
cates that the lowest Matsubara mode contains 
the non-perturbative physics js^,^. It therefore 
makes sense to perturbatively integrate out the 
higher Matsubara modes, and obtain a 3D effec- 
tive action for the static mode. This can then be 
treated by lattice methods. 

The final result for the pressure takes the form 
of an integral over a coupling constant in the 3D 
theory, with a parameter cq associated with the 
constant of integration which can be determined 
by an arduous 4-loop calculation in the 4D theory. 

In Fig. 1^ we see the range of predictions for the 
pressure of hot QCD with no quarks for various 
values of the constant cq. The value eo — 10(1) 
agrees best with 4D lattice calculations. The ef- 
fort to calculate cq independently is currently un- 
derway. 

This method is certainly better behaved than 
4D perturbation theory. It may also avoid some of 
the difficulties that beset 4D lattice calculations. 
For example, light quarks only come in through 
the coefficients in the 3D action, and so do not im- 
pose any additional computational burden. Also 
it is possible to perform calculations at /i > . 
But it also inherits some of their limitations, such 
as being restricted to thermodynamic observables 
that can be measured in Euclidean time. 
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Figure 4. The pressure of hot gluons, calculated 
by dimensional reduction, for various values of the 
undetermined constant eo [Q. Statistical errors 
are shown for e = 10. The dotted line indicates 
predictions of lattice calculations. 



5. Lattice QCD at high temperature and 
low density 

Lattice calculations at non-zero density (chem- 
ical potential > 0) have to contend with the 
"sign problem". Until recently this meant that 
there were no lattice results for realistic (3 color, 
3+1 dimensional) QCD away from the /z = axis 
of the phase diagram. There has been a lot of 
work on more tractable variants such as QCD 
with two colors ||8^, QCD at finite isospin den- 
sity Gross-Neveu models in 2-1-1 dimensions 
pTI] and so on. It has also been suggested that at 
asymptotically high density the sign problem in 
QCD might become more tractable js^ . 

In the last few years this situation has changed, 
with three groups reporting lattice results for the 
position of the crossover line and even the criti- 
cal point in the high-T low-/i region of the QCD 
phase diagram. Before discussing those results, I 
will review the nature of the sign problem. 



The sign problem is not a deep problem of 
principle, but it is a very inconvenient feature 
of the standard formulation of the QCD func- 
tional integral. The problem is that when ^ 7^ 0, 
the factor of det(Af[A]) in the gauge field func- 
tional integral is no longer a positive real 
number. This means that the only practical way 
of evaluating such a high-dimensional integral, 
the Monte-Carlo importance sampling method, 
which interprets the Boltzmann weight B[A] = 
det(M[yl]) exp(— S'gauge) as a probabihty, cannot 
be applied. 

In principle it might seem that the sign problem 
could be straightforwardly fixed by "reweight- 
ing" One example of this would be to split 

off the phase of the weight by writing B[A] = 
\B[A] \ exp{i(j)[A]) , and use the absolute value of 
the Boltzmann factor for importance sampling. 



jm\B[A]\ 



(0)|| = {l/Z^\)JmO[A]\B[A]\ 



(6) 



The true expectation value of the observable can 
be obtained by including the complex phase in the 
observable, taking its expectation value in the \B\ 
ensemble, and dividing by a reweighting factor 



{0) = {OeM^m))\\/iZ/Zu) 



(7) 



But the reweighting factor is the ratio of two par- 
tition functions, so it is related to the difference 
between their free energy densities / 



-exp(~V^(/~/||)). 



(8) 



The true partition function Z contains the phase 
factor exp(i(/)[^]), and is therefore smaller than 
Z\, so / > /||. Thus the reweighting factor in 



|) becomes exponentially small in the thermody- 
namic limit V — > 00. The same is true for (0)|[, 
because (O) itself is not exponentially large in V. 
This means exponentially large statistics are re- 
quired to obtain any given accuracy in (O). 

It must be stressed that the sign problem arises 
from the way we structure our evaluation of the 
functional integral, integrating out the fermion 
fields for each fixed gauge field configuration. The 
physically relevant quantities such as Z and (O) 
are real and positive, and if we could restruc- 
ture the ensemble of gauge field configurations 
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into families within which the weights and ob- 
servables analytically added up to positive values 
then there would be no sign problem. This is 
what cluster algorithms can do, and such a re- 
structuring has been achieved for the Z3 model 
that approximates QCD at very high quark mass 
and chemical potential |Q . One strategy for the 
ultimate solution of the sign problem is to recast 
the QCD functional integral into a form to which 
cluster algorithms can be applied . 

In the meantime, it makes sense to try to per- 
form calculations using the methods available, in 
regions of the phase diagram where the sign prob- 
lem is weakest. This is the region of small /x and 
high temperature. In the last few years, various 
groups have have been exploring techniques ap- 
plicable at /X < T: multi-parameter reweighting, 
derivatives with respect to at /i = 0, and imag- 
inary chemical potential. 
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Figure 5. The position in the /i-T plane of 
the crossover between the hadronic and QGP 
phases, including an estimate of the position of 
the critical point at {fi*,T*) and the first-order 
line, obtained from a lattice calculation with two- 
parameter reweighting Note that the hori- 
zontal axis is /is = 3/x. 



5.1. Multi-parameter reweighting. 

Fodor and Katz |Q suggested a reweighting 
algorithm like the one described above, but using 



the /i = ensemble to measure observables rather 
than the \B\ ensemble. They exploit the freedom 
to vary the temperature of the = ensemble, 
and thereby improve its overlap with the desired 
^ > ensemble. This only postpones the sign 
problem to larger volumes, but from their initial 
results it appears that reasonable results may now 
be obtained on lattices that are large enough to 
be physically interesting. Their estimate of the 
crossover line and chiral critical point is shown in 
Fig. H (for the equation of state see Q). The nu- 
merical values should be interpreted with caution: 
the light quark masses are well above their phys- 
ical values, so the true critical point is presum- 
ably at a somewhat lower chemical potential . 
The lattice quark action is unimproved (Kogut- 
Susskind) and used at a relatively coarse lattice 
spacing a ^ 0.2 fm where it is known to have large 
discretization errors. Thus there are potentially 
large systematic errors, and as with all reweight- 
ing methods the statistical errors can be mislead- 
ingly small. Nonetheless, this is the only method 
that has obtained the position of the chiral criti- 
cal point, and it seems there is a real chance that 
it may be able to do so for physical quark masses. 
This is a striking improvement over the situation 
two years ago. 

5.2. d/dfi at /X = 0. 

A closely related strategy for extending lattice 
QCD to small chemical potential is to evaluate 
derivatives of observables with respect to /x at 
/X = 0, where there is no sign problem iQ. This 
is like taking the /x-derivative of Fodor and Katz's 
reweighting factor, and including that in the ob- 
servable. The relevant calculations have been 
performed by the Bielefeld/ Swansea group [ p7[ . 
Their observables are the susceptibilities for the 
Polyakov loop and chiral condensate, the max- 
ima of which indicate the position of the crossover 
line in the /i-T plane. Using improved gauge and 
staggered quark actions on a lattice of spacing 
a « 0.2 fm, they obtain the results shown in 
Fig. ^. The results agree very well with the multi- 
parameter reweighting lattice calculation, whose 
estimate of the critical endpoint is also shown. 
It remains to be seen whether this approach has 
advantages over multi-parameter reweighting. In 
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particular, it is very hard to know how deep into 
the ^ > region the Taylor expansion around 
/i = will continue to work. To the left of the 
dotted n/T = 0.4 hue in Fig. || the sign problem is 
under control, but the importance of higher order 
/i-derivatives is yet to be checked. 
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Figure 6. QCD hadronic/QGP crossover line 
from derivatives of susceptibilities with respect 
to ^, evaluated at = Note the excellent 
agreement with Fodor and Katz's multiparame- 
ter reweighting calculation of the critical point 
(/i*,T*) (diamond). 



5.3. Imaginary /i. 

Over the years, various workers have pointed 
out that the sign problem only occurs for real 
chemical potential — the functional integral with 
imaginary /i is amenable to standard Monte-Carlo 
methods . Now full 4D calculations by de For- 
crand and Philipsen |Q have yielded estimates of 
the position of the critical line in the /i-T plane 
for two-flavor QCD (for early explorations of 4 
flavor QCD, see |10C| ). These authors perform 
calculations at several imaginary values of /i, and 



then analytically continue to real /i by fitting their 
results to a power series in and extrapolat- 
ing from negative to positive ^j? [101|. The ob- 



servables used to identify the crossover line were, 
as usual, the susceptibilities for the plaquette, 
Polyakov loop, and chiral condensate. Above the 
critical point (/i > fX*) there is a line in the jji-T 
plane along which these diverge in the large vol- 
ume limit: this is the first-order line. Below the 
critical point (/i < /i*) each susceptibility shows a 
"ridge" of finite height in the large volume limit, 
and these ridges are expected to lie close together, 
approximately defining a narrow crossover strip. 
At large but finite volume, the line defined by 
the ridge/divergence of each susceptibility x "^iH 
be an analytic function T^{p). At infinite vol- 
ume different susceptibilities will give the same 
line above fj,* but may give different ones below 
fi* so it is not clear that the function remains an- 
alytic for all /i, but at the volumes studied so far 
r^(/i) turns out to be a very docile expansion in 
/i^ for imaginary /j,, fitting very well to the form 
Co + ci/i^ with no need for a /i^ term. The resul- 
tant crossover line for QCD with two light flavors 
is plotted in Fig. |^. It is consistent with the re- 
sults from the other methods, and the errors are 
smaller. 

Imaginary chemical potential is quite different 
from the two methods described above, and the 
errors appear to be well under control. It will 
be very interesting to see whether it can estimate 
the location of the critical point as accurately as 
it delineates the crossover region. 

5.4. Relevance to experiment 

It is very impressive that the different methods 
described above, all of which are being applied 
to QCD for the first time, achieve agreement in 
the position of the crossover line, and that there 
is even a rough estimate of the position of the 
critical point fi*g - 700 MeV, T* - 160 MeV. 
Putting aside caution and taking the numbers se- 
riously for a moment, could such a critical point 
be accessible to heavy ion collision experiments? 
Treating the estimated value of as an upper 
limit, and bearing in mind that a heavy ion col- 
lision scans a range of temperatures down to the 
chemical freeze-out temperature Tfo as the fire- 
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Figure 7. QCD hadronic/QGP crossover line 
from analytic continuation of imaginary fi re- 
sults to real /i |9^. The result is consistent 
with Figs. Note that the horizontal axis is 

flB = 3fJ,. 



ball expands, we can see from Table |l| that there 
is a chance that quark matter near the critical 
point could have been made at SPS. It may also 
be produced at the future experiment planned at 
GSI [102|. And now we can even envisage the pos- 
sibility that theorists may predict the position of 
the critical point before experimentalists observe 
it. 



experiment 




Tfo 


RHIC 


40 MeV 


170 MeV 


SPS 


250 MeV 


170 MeV 


AGS 


550 MeV 


130 MeV 


SIS 


700 MeV 


80 MeV 



Table 1 

Approximate chemical freeze-out temperatures 
and baryon chemical potentials for heavy-ion col- 
lisions at various colliders [103|. 



high temperature is their application to heavy 
ion collisions at current and future heavy-ion col- 
liders. The success of hydrodynamic models of 
RHIC collisions [ 104 indicates that equilibrium 
is achieved within a very short time ^ 1 fm/c, 
but perturbative estimates give longer timescales 
\10l\. It is therefore very important to under- 
stand nonperturbative field theory far from equi- 
librium. 

In the last two years there has been signifi- 
cant progress in treating non-equilibrium scalar 
field theory. It is based on systematic loop 
1 10^ and 1/iV |10|,|109| expansions of the two- 
particle-irreducible (2PI) effective action formal- 
ism p^ , |65[ discussed in section 4.2. The formal- 
ism is illustrated diagramatically in Fig. |[ A 
typical mean-field calculation would only include 
the 0{N) tadpole (momentum-independent) con- 
tribution to the self-energy, and would not show 
thermal equilibration at long times. The NLO 
2PI calculation also includes the 0(1) bubble 
chain contribution of Fig. |^, which is crucial for 
including rescattering effects that lead to equili- 
bration. 




> 0(1) 



6. Far from equilibrium 

One motivation for calculations of the QCD 
equation of state and transport properties at 



Figure 8. Coupled NLO Schwinger-Dyson equa- 
tions for the full propagator from the 2PI effective 



action approach of [107 
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Fig. H shows the time evolution of the propa- 
gator G{p) in 1+lD scalar field theory, for three 
momentum modes |p| = 0,3,5. The solid lines 
are for initial conditions close to the equilibrium 
value, the dashed lines for initial conditions fur- 
ther from the equilibrium value, and the dot- 
dashed lines are for initial conditions far from the 
equilibrium value. 



The rate of approach to equilibrium 
is found to be quite different from the initial 
damping rate 'y('^™p), and the overall equilibra- 
tion time in scalar theories is clearly much longer 
than the microscopic timescales. Obviously it 
will be very interesting to push these calcula- 
tions to higher dimensions [110| where sponta- 
neous symmetry breaking can be studied, and 
to extend these methods to gauge theories, to 
see whether QCD predicts equilibration times as 
short as 1 fm/c for the QGP. 
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Figure 9. Time dependence of three modes of 
the equal time propagator in 1-l-lD 0* theory at 
non-zero temperature, for various initial condi- 
tions. Calculated by 3- loop 2PI effective action 
1 108 1 . Note the quickly damped initial oscilla- 



tions, followed by drifting and then comparably 
slow approach to equilibrium. 



We see that the modes evolve in three stages, 

1. Early: Initial oscillations quickly damped 
with rate j('^'"^p\ 

2. Intermediate: Drifting towards the equilib- 
rium value with non-exponential/power-law 
behavior. 

3. Late: Exponential approach to equilibrium, 
with rate 7(tii°™) ^ 7(damp)_ 



7. Conclusions 

There is great enthusiasm and activity among 
theorists studying high-temperature and high- 
density QCD, stimulated by the successes of cur- 
rent heavy-ion experiments and the anticipation 
of future ones, and also by ongoing advances 
in the quality of observational data for neutron 
stars. 

In high-density QCD, the upsurge of interest in 
color superconductivity is maturing into a effort 
to develop signatures for the presence of color- 
superconducting quark matter in compact stars. 
While ultra-high-density first-principles calcula- 
tions and medium-density NJL calculations are 
becoming more sophisticated, useful results also 
follow from parameterizing the quark pairing with 
a gap parameter and studying its effects on phe- 
nomenologically important quantities. 

In high-temperature QCD, lattice gauge the- 
ory is making continuous advances, achieving re- 
liable calculations of equilibrium properties of the 
QGP such as screening lengths and the equation 
of state at temperatures up to a few times Tc- Di- 
mensionally reduced lattice calculations will be 
able to extend this to arbitrarily high tempera- 
tures, when one remaining constant is computed. 
Various resummed perturbative calculations ap- 
pear to agree with the lattice on the approximate 
size of deviations from the free equation of state 
above the critical temperature, although they are 
not yet on a solid independent footing. Complete 
leading-order perturbative expansions of trans- 
port properties, including all leading logs, are be- 
ing achieved. The photon emission rate of the 
QGP has been computed, and other quantities 
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wi ll soon be calculated. 

jLattice QCD has at last taken its first step 



E. lancu, P hys. Rcpt. 359 (2002) 355 jiep 



into the non-zero-density region of the phase di- 
agram. Remarkably, there is already agreement 
among the different methods on the position of 
the crossover line, and there is even a first esti- 
mate of the position of the chiral critical point, 
which is a sharply-defined physical observable. It 
would therefore be valuable for theorists, as well 
as striving for more accurate predictions of the 
position of the critical point, to further develop 
the existing work on signatures of the critical 
point in heavy-ion collisions. 

Finally, field theories far from equilibrium 
are becoming accessible to controlled calcula- 
tion. The 2PI effective action method has ap- 
plications to disoriented chiral condensates, crit- 
ical fluctuations, and other complexities of heavy 
ion collisions, as well as cosmological questions 
such as the exit from inflation in the early uni- 
verse. Efforts to apply this method to higher- 
dimensional scalar field theories, spontaneous 
symmetry breaking, and gauge theories are in 
progress. 

Strong interaction physics at high temperature 
and density is a promised land that lies at the 
heart of modern physics, bordering on particle 
physics, nuclear physics, astrophysics, and con- 
densed matter physics. Our first explorations are 
revealing the true extent of its fertility and rich- 
ness. 
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